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Abstract 

We report the first observation of a b — > u type charmless baryonic B decay, B + — ► ppir + , as 
well as 6 — > s type i?° — ► ppK° and -B + — > ppK* + decays. The analysis is based on a 78 fb _1 
data sample recorded on the T(4S) resonance with the Belle detector at KEKB. We find B(B + — ► 
PPTT+) = (3.0618;^ ± 0.37) x KT 6 , B{B° -► ppK°) = (1.88±g;^ ± 0.23) x KT 6 , and B{B+ -► 
ppi^*+) = (10.3±|;|±i;?) x lO" 6 . We also update B(B+ -► ppA" + ) = (5.66t^ ± 0.62) x 10~ 6 , and 
present an upper limit on B(B° — > ppK*°) at the 90% confidence level. A common feature of the 
observed decay modes is threshold peaking in baryon pair invariant mass. 

PACS numbers: 13.25.Hw, 13.60.Rj 
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The Belle Collaboration recently reported the observation of B + — > ppK + 1] and B° — > 
decays. These are the first examples of B meson decays to charmless three-body 
final states containing baryons, and are candidates for b — > s penguin transitions. Our 
observation of these modes has stimulated much theoretical interest |E HI IE IE IE El 
due to these decay channels may be used to test our theoretical understanding of rare decay 
processes involving baryons and search for direct CP violation. One interesting feature of the 
observation is that the baryon pair mass spectra seem to peak toward threshold as originally 
conjectured by Refs. [IE U|- Lately more speculations on this are proposed 0, IE IE] • m 



this letter we report the first observation of B + — ► ppn + which is dominated by the 



b — > u tree diagram. We also report the first observation of B° — > ppK° and B + — >• ppK* + 
decays, and improve the measurement of B + — ► ppK + . A search for the B° — > ppK*° mode 
yields only an upper limit. With these new results we study the pp mass spectra to see if 
the threshold peaking observed in our previous results is confirmed with the newly observed 
modes. 

We use a 78 fb _1 data sample, consisting of 85.0 ± 0.5 million BB pairs, collected by the 
Belle detector at the KEKB asymmetric energy e + e~ (3.5 on 8 GeV) collider |13| . The Belle 
detector is a large solid angle magnetic spectrometer that consists of a three layer silicon 
vertex detector (SVD), a 50 layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time of flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of Csl (Tl) crystals located inside 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return 
located outside of the coil is instrumented to detect K\ mesons and to identify muons. The 
detector is described in detail elsewhere 3- 



The event selection criteria are based on the information obtained from the tracking 
system (SVD+CDC) and the hadron identification system (CDC+ACC+TOF), and are 
optimized using Monte Carlo (MC) simulated event samples. All primary charged tracks are 
required to satisfy track quality criteria based on the track impact parameters relative to the 
interaction point (IP). The deviations from the IP position are required to be within ±1 cm in 
the transverse (x-y) plane, and within ±3 cm in the z direction, where the z axis is defined 
by the positron beam line. Proton, kaon and pion candidates are selected using p/K/n 
likelihood functions obtained from the hadron identification system. For protons, we require 
L p /(L p + L K ) > 0.6 and L p /(L p + L n ) > 0.6, where Lp/x/^ stands for the proton/kaon/pion 
likelihood. We require L K j[L K + L n ) > 0.6 to identify kaons and L n /(L K + L n ) > 0.6 for 
pions. Kg candidates are reconstructed via the tt + tt~ decay channel and have an invariant 
mass with \M n + n - — M K o\ < 30 MeV/c 2 . The candidate must have a displaced vertex 
and flight direction consistent with a K s originating from the interaction point. We use the 
selected kaons and pions to form K* + (— > K s tt + ) and K*° (— > K + it~) candidates by requiring 
the invariant mass of the combination to be within 80 MeV/c 2 of the K* mass. The numbers 
used are 892 MeV/c 2 and 896 MeV/c 2 for M K *+ and M K *o, respectively. To ensure the decay 
process is genuinely charmless, we apply the charm veto. The regions 2.85 < M p p < 3.128 
GeV/c 2 and 3.315 < M pP < 3.735 GeV/c 2 are excluded to remove background from modes 
with r] c , J /if) and if)', Xco, Xd mesons, respectively. The region 2.262 < M pK o < 2.310 GeV/c 2 
is excluded to remove A c candidates. 

Since the initial collision center of mass energy is set to match the T(4S) resonance, 
which decays into a BB pair, one can use the following two kinematic variables to identify 



the reconstructed B meson candidates: the beam constrained mass Mb c = \J E^ eam — p 2 B , 
and the energy difference AE = E B — E hcam , where i?beam, Pb, and E B are the beam energy, 
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the momentum, and energy of the reconstructed B meson, respectively, in the rest frame 
of the T(4S). The candidate region is defined as 5.2 GeV/c 2 < M bc < 5.29 GeV/c 2 and 
\AE\ < 0.2 GeV. 

The dominant background arises from the continuum e + e~ — > qq process, with much 
smaller contributions from "cross-feed" , where similar types of rare decay events pass each 
other's signal criteria. The background from charm-bearing and charmless mesonic decays 
is negligible. In the T(4S) rest frame, continuum events are jet-like while BB events are 
more spherical. One can use the reconstructed momenta of final state particles to form 
various shape variables (e.g. thrust angle, Fox- Wolfram moments, etc.) in order to catego- 
rize each event. We follow the scheme defined in Ref. ^| that combines seven event shape 
variables into a Fisher discriminant 0] in order to suppress continuum background. The 
variables chosen have almost no correlation with M bc and AE. Probability density func- 
tions (PDFs) for the Fisher discriminant and the cosine of the angle between the B flight 
direction and the beam direction in the T(4S) rest frame are combined to form the signal 
(background) likelihood £s(bg)- The signal PDFs are determined using signal MC simula- 
tion; the background PDFs are obtained from the continuum MC simulation for events with 
5.27 GeV/c 2 < M bc <5.29 GeV/c 2 and \AE\ < 0.1 GeV. We require the likelihood ratio 
CTZ = C S /(C S + C BG ) to be greater than 0.8, 0.9, 0.85, 0.8 and 0.8 for ppK + , pp7r + , ppK° s , 
ppK* + and ppK*° modes, respectively. The selection points are determined by optimization 
of S/ (S + B), where S and B denote the number of signal and background, respectively. 
Note that a nominal signal branching fraction of 4 x 10~ 6 is assumed. 

Since the previous studies 0, 0| showed an enhancement at low baryon-antibaryon in- 
variant mass, we first focus on the near threshold region by requiring M p p < 2.85 GeV/c 2 
to make sure it is below charmonium threshold. The M bc distributions (with \AE\ < 0.05 
GeV), and the AE distributions (with M bc > 5.27 GeV/c 2 ) for the ppK + , ppTr + , ppKg, and 
ppK* + modes are shown in Fig. d We use an unbinned likelihood fit to estimate the signal 
yield: 

N 

L = Y[[f s P s {M hCi , AEi) + (1 - f s )P b (M hc ., AEi)), 

where P s {Pb) denotes the signal (background) PDF and f s is the signal fraction of the total 
iV candidates. For the signal PDF, we use a Gaussian in M bc and a double Gaussian in 
AE. We fix the parameters of these functions to values determined by MC simulation. 
Background shapes are studied using sideband events: 0.1 GeV < \AE\ < 0.2 GeV for the 
M bc study and 5.20 GeV/c 2 < M bc < 5.26 GeV/c 2 for AE. These shapes are confirmed with 
a continuum MC sample. We use the following parametrization first used by the ARGUS 
collaboration, /(M bc ) oc M bc \/1 ~ % 2 exp[— £(1 — x 2 )}, to model the M bc background, where 
x is defined as M bc /.E beam and £ is a parameter to be fit. The AE background shape is 
modeled by a first order polynomial. There are possible cross-feeds from ppK* + and ppK*° 
modes to ppK + , ppii + and ppK^ modes; therefore the cross-feed region (AE < —0.16 GeV) 
is excluded in the fit. Since the ppn + mode can contain non-negligible cross-feed events from 
the ppK + mode, we include the ppK + MC cross-feed shape in the fit for the determination 
of the ppir + yield. The fit results are shown in Fig. ^ by solid curves. The fit yields are 96.4 
t\H, 37 A ±$, 11.3 til, and 14.5 tto with significances of 15.3, 6.7, 5.1, and 6.0 standard 
deviations for the ppK + , ppn + , ppK®, and ppK* + modes, respectively. The significance is 

defined as J— 2\n(L / L max ) [13], where L and L max denote the likelihood with signal yield 
fixed at zero and at the fitted value, respectively. 
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Although we do not have adequate statistics to perform a full Dalitz plot analysis, the 
observed distribution is not uniform over phase space. To reduce the model dependence in 
determining the branching fraction, we fit the signal yields separately in nine bins of M pp , 
and correct for the detection efficiencies from MC simulation in each bin. In Fig. |2J we show 
the signal yield versus M pp , with three-body phase space from MC (normalized in area) 
superimposed. The observed mass distributions all peak at low pp mass. The branching 
fractions (B) in bins of M pp for the observed modes are given in Table The upper limit of 
the last bin is different for each mode and is equal to the kinematic limit, M P p_n m . We sum 
the partial branching fractions to obtain the total branching fractions. The results are listed 
in Table |n]and the branching fractions below charmonium threshold, M pp < 2.85 GeV/c 2 , 
are also listed for comparison. Note that B(B° — > ppK°) = 2B(B° — > ppKg) is assumed. 

The search for ppK*° gives a yield of 13_5 events with a significance of about 3 standard 
deviations. Since it is less significant, we use the fit results to estimate the expected back- 
ground and compare this with the observed number of events in the signal region in order to 
set the upper limit on the yield at the 90% confidence level 3, 111- Note that the systematic 



uncertainty is needed in this estimation. The upper limit yield of the ppK*° mode in the 
full M pp range is determined to be 57 at the 90% confidence level. The branching fraction 
is found to be B(B° ppK*°) < 7.6 x 10" 6 . 

Systematic uncertainties are studied using high statistics control samples. For proton 
identification, we use aA^ pn~ sample, while for K/ir identification we use a D* + — > D°7T + , 
D° — > K~n + sample. Tracking efficiency is studied with fully and partially reconstructed 
D* samples. Kg reconstruction efficiency is studied with a D~ — > K s 7t~ sample. The 
£71 continuum suppression uncertainty is studied with B + — > J/ipK + , B° — > J/ipK s , 
B + — > J/ipK* + , and B° — > J/ipK*° (with J/ip — > control samples. Based on these 

studies, we assign a 1% error for each track, 3% for each proton identification, 2% for each 
kaon/pion identification, 5% for K s reconstruction and 6% for the CR, selection. 

The systematic uncertainty in the fit yield is studied by varying the parameters of the 
signal and background PDFs. We assign an error of 5% for the ppn + mode and 4% for the 
other modes. There is a possibility of non-resonant Kgir + combinations passing our K* + 



TABLE I: Branching fractions (in units of 10 6 ) in different M pp bins for the ppK + , pp?r + , ppKg 
and ppK* + modes. 



M pp (GeV/c 2 ) 


ppK + 


ppiT + 


ppK° s 


ppK* + 


1.876 - 2.0 


91+ - 23 

u - yl -0.20 


29+ - 15 


o.38±g:iS 


i ' 4U -0.74 


2.0-2.2 


i OO+0.31 
l.OO_ 28 


0.57t°;i 


0-17^ 


2 11+ 1 ' 17 


2.2 - 2.4 


1 20+ - 29 


0.39j&g 


o.oo^i 2 


i m +1 - 00 


2.4-2.6 


0.82±°;| 


0.461^ 


o.oe^o 2 


64+ - 91 


2.6-2.8 


^i+ - 21 

u - oi -0.17 


o.os^i 2 


O-OS^o 2 


Q^+ - 94 
u - yo -0.60 


2.8-3.4 


^+ 22 
u - o,5 -0.18 


o oo +0 - 04 

u - uu -0.04 


o u +0 - 15 


1.19i ;?? 


3.4-4.0 


1^+ 015 

u - i0 -o.n 


u - io -0.11 


o oo +0 - 07 

u - uu -0.07 


2.0711;! 


4.0-4.6 


o on +008 

u - uu -0.08 


u - io -0.21 


0-io±8;i8 


0-93±i;?l 


4.6 - M p p_i im 


20+ 17 

u - zu -0.12 


J-.UU_q.45 


o.oo^il 


n/a 
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TABLE II: Branching fractions (in units of 10 6 ) in the full M pp range and below the charm 
threshold (M pp < 2.85GeV/c 2 ) for the ppK+, ppir + , ppKg and ppK* + modes. Statistical and 
systematic errors are quoted. 





mode full M pp range 


M pp < 2.85GeV/c 2 


B + 


-^ppK+ 5.66±g;|? ±0.62 


4.89±gJ§ ± 0.54 


B+ 


-^ppir+ 3.06±g;^ ±0.37 


1.76±g;^ ± 0.21 


B° 


-^ppK° 1.88^;^ ±0.23 


1.56lg;g ± 0.19 


B+ 


-» ppK*+ 10.3ii|+i.3 


fi 7+2.4+0.9 
u -' -2.0-1.1 



cut. We study this by fitting for the B yield in bins of M K o^+. Choosing a Breit-Wigner 
signal and a first order polynomial background to estimate the B + — > ppK* + branching 
fraction, the number drops by 9%. We include this uncertainty in the estimation of signal 
yield for the ppK* + mode. The MC statistical uncertainty and modeling with nine M pp bins 
contributes a 2% error in the branching fraction determination. The error on the number of 
BB pairs is determined to be 1%, where an assumption is made that the branching fractions 
of T(4S) to neutral and charged BB pairs are equal. 

We first sum the correlated errors linearly (e.g., total 5% tracking error for the ppK* + 
mode) and then combine with the uncorrelated ones in quadrature. The determined system- 
atic uncertainties are 11%, 12%, 12%, and 13% for the ppK + , ppn + , ppK°, ppK* + , 
and ppK*° modes, respectively. 

Using the charm veto events, we perform a cross-check of our analysis. We deter- 
mine B — > J/ipK^ branching fractions by following the same analysis procedure with 
3.07 GeV/c 2 < M pp < 3.11 GeV/c 2 and using B(J/i[) -> pp) = (2.12 ± 0.10 x 10~ 3 ) 
0. The obtained branching fractions are B{B + -> J/ipK + ) = (1.17±g^| ± 0.16) x 10" 3 , 
B(B° -> J/^jK°) = (1. 16^ ±0.16) x lO" 3 , B{B + -> J/^K*+) = (1.08±$;| ± 0.17) x lO" 3 , 
and B{B° -* J/i/jK* ) = (i.40±o;ie ± °- 21 ) x 10_3 ' which are in agreement with the world 
average values [171 ]. 

The present results, shown in Table [HI offer valuable information for understanding the 
mechanism of charmless baryonic B decay. In particular, the threshold peaking behavior 
is now firmly established. The B + — ► ppK + data can be used to constrain the production 
of narrow glueball states that decay to pp 0. With the current statistics, we can not set 
a very stringent bound. It should be noted that the observed B + — > ppir + rate is less 
than B + — > ppK + , which is consistent with what is observed in B — > Ktt, tttt modes. 
The B° —>■ ppK° rate is considerably lower than that of the B + — > ppK + mode which 
should be contrasted with B°' + — > tt°K 0,+ modes and also B 0,+ — > J/i/jK 0,+ modes . This 
indicates that the intermediate pp system is nontrivial. These modes are of interest for direct 
CP violation searches. For the B^ — ► ppK T and pp7t T modes that have larger statistics, 
we define the charge asymmetry as (N B - — N B +)/(N B - ± N B +) and find the values are 
—0.05 ±0.11 ±0.01 and —0.16 ±0.22 ±0.01 respectively. The systematic error is determined 
by checking the null asymmetry with a B — > D(— > K 71)71 control sample. The measured 
asymmetries are consistent with zero for their large statistical uncertainties. 
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FIG. 1: M bc and AE distributions for (a) ppK + , (b) ppir + , (c) ppK$, and (d) ppK* + modes, for 
M pp < 2.85 GeV/c 2 . 
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FIG. 2: Fitted signal yield divided by the bin size for (a) ppK + , (b) ppir + , (c) ppK®, and (d) 
ppK* + modes in bins of M p p. The shaded distribution is from the phase-space MC simulation with 
area normalized to signal yield. 
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